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1 SCOPE 


The  objective  of  this  paper  is  to  describe  as  completely  as  is 
practical  the  present  state-of-the-art  for  predicting  the  structural, 
response  of  surface-loaded  pavements.  In  scope  the  paper  primarily 
focuses  on  the  finite  element  techniques  as  related  to  pavement  analysis. 
Other  techniques  and  topics  are  introduced  to  provide  a complete  picture 
and  offer  a context  for  discussion. 

Most  of  the  information  presented  is  based  on  the  personal  ex- 
periences of  the  authors  and  their  colleagues.  No  attempt  has  been  made 
to  compile  an  exhaustive  list  of  computer  codes  applicable  to  pavement 
analysis  and  design.*  However,  each  area  relating  to  state-of-the-art 
prediction  of  pavement  response  is  discussed  and  summarized. 

Common  to  all  prediction  techniques  are  the  fundamental  laws  of 
mechanics  as  typified  by  the  theory  of  elasticity.  Therefore,  to  this 
extent  all  prediction  techniques  share  a common  origin.  However,  this 
commonality  becomes  obscured  as  soon  as  the  assumed  mathematical  model 
is  defined.  That  is,  by  specifying  geometry,  constitutive  laws,  and 
boundary  conditions,  a plethora  of  prediction  techniques  emerge  ranging 
from  simple  linear  plate  theory  to  complex,  nonlinear,  three-dimensional 
continuum  theory.  Moreover,  several  techniques  applied  to  a common 
situation  may  produce  different  results.  This  leads  to  confusion  and 
provides  a challenge  for  this  paper  to  meet. 


• A list  of  pavement  analysis  techniques  has  been  compiled  by  Laboratoire 
Central  des  Ponts  et  Chaussees.  These  techniques  are  mostly  elastic 
layer  and  finite  element.  Most  of  the  techniques  listed  are  in  use 
outside  the  United  States. 
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2 THEORETICAL  BASIS  OF  THE  PRINCIPAL 
ANALYTICAL  TECHNIQUES 

The  techniques  considered  are  Westergaard,  elastic  layer  (i.e., 
the  Burmister  problem,  including  Boussinesq),  and  finite  element.  These 
are  the  principal  techniques  used  for  pavement  analysis.  Each  of  these 
methods  is  based  on  the  fundamental  lavs  of  mechanics  which  require  con- 
sideration of  four  distinct  mechanistic  concepts:  (a)  equilibrium 

equations,  (b)  kinematical  relations,  (c)  constitutive  lavs,  and  (d) 
boundary  conditions. 

2.1  MATHEMATICAL  CONCEPTS  FOR  ALL 
TECHNIQUES 

All  techniques  require  that  the  "real"  pavement  be  idealised  by 
a mathematical  model.  The  use  of  mathematics  requires  definition  of  a 
domain  to  be  simulated  by  the  idealisation.  The  "principal"  techniques 
present  two  choices:  the  semi-infinite  domain  characteristic  of 

Westergaard  and  elastic  layer,  and  the  finite  domain  or  "soil  island" 
typical  of  finite  element.  Usually  pavement  problems  are  most  easily 
modeled  as  semi-infinite  bodies,  and  the  necessity  for  truncation  of  the 
"real"  situation  is  an  obstacle  in  the  application  of  the  finite  element 
technique. 

2.1.1  EQUILIBRIUM 

Static  equilibrium  equations  thematic ally  express  the  concept 

that  the  summation  of  forces  at  any  point  vlthin  the  idealisation  and 
on  the  boundaries  of  the  idealization  must  equal  zero.  For  the  classical 
techniques  used  to  solve  the  Westergaard  and  elastic  layer  idealizations, 
this  requires  the  solution  of  partial  differential  equations.  These 
equations  are  usually  complex  enough  that  they  must  be  solved  by  computer. 
While  the  equations  are  an  "exact"  expression  for  equilibrium,  use  of 
numerical  procedures  to  solve  them  does  introduce  some  level  of  approxi- 
mation. On  the  other  hand,  equilibrium  equations  derived  by  using  the 
finite  element  method  are  a linear  algebraic  set  of  equations  that 
approximate  the  equilibrium  condition  for  the  structure.  While  the 
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number  of  equations  and  coefficients  is  large,  they  are  readily  formed 
and  solved  by  the  computer. 

While  finite  element  techniques  are  generally  not  as  exact  as 
classical  ones,  their  versatility  and  ease  of  derivation  are  vastly 
superior  in  most  cases.  For  a limited  range  of  problems  (e.g.,  linear, 
layered  pavements),  classical  solutions  offer  the  most  efficient  way  to 
calculate  responses  at  a few  points;  however,  for  a large  number  of 
points,  finite  element  techniques  are  much  more  efficient.*  This  follows 
from  the  requirement  that  for  a classical  idealization,  the  partial 
differential  equations  must  be  resolved  at  each  location  where  responses 
are  computed.  Conversely,  the  finite  element  technique  solves  for  all 
responses  of  the  idealization  at  once. 

In  summation,  study  of  the  equilibrium  concept  shows  that  clas- 
sical solutions  in  limited  circumstances  may  be  more  computationally 
efficient  than  finite  element  schemes;  however,  a high  price  in  versa- 
tility is  paid  for  the  use  of  clrssieal  solutions. 

2.1.2  KINEMATICAL  RELATIONS 

Another  part  of  any  mathematical  idealization  of  a pavement  is 
the  relation  of  displacement  to  strain  and  the  variation  of  ctrain- 
displacement  within  the  idealization.  A linear  strain-displacement 
relation  for  most  airfield  idealizations  is  appropriate  (i.e.,  small 
displacement  theory  may  be  used).  An  exception  to  this  occurs  when 
significant  deflections  and  membrane  forces  are  present  in  a pavement 
layer;  an  example  is  shown  in  Section  L.5.  A more  significant  aspect 
of  kinematical  relations  is  the  consideration  of  discontinuous  strain- 
displacement  fields  to  model  the  phenomena  of  interlayer  debonding  and 
Joints. 

Implementation  of  complicated  kinematical  relations  highlights 
the  distinctive  advantage  of  numerical  procedures.  While  some  compli- 
cation in  .we  numerical  algorithm  results,  no  significant  increase  in 
difficulty  for  the  solution  of  the  equilibrium  equations  occurs.  Intro- 
duction of  nonlinear  discontinuous  displacement  fields  in  the  partial 
differential  equations  of  the  classical  techniques  in  all  but  the 
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simplest  circumstances  produces  an  intractable  problem.  In  summary, 
for  kinematically  complicated  problems  only  numerical  procedures  such 
as  finite  element  are  of  practical  benefit. 

2.1.3  CONSTITUTIVE  LAWS 

For  the  prediction  of  pavement  response,  specification  of  appro- 
priate constitutive  laws  is  the  only  significant  challenge  still  to  be 
solved  conceptually.  Constitutive  equations  mathematically  represent 
the  observed  phenomenological  relations  between  stress  and  strain  in  a 
given  material  and  are,  at  best,  approximations  of  the  'Veal"  material 
behavior.  Materials  axe  characterized  using  a variety  of  measures  (for 
example,  bulk,  shear,  and  Young's  modulus;  Poisson's  ratio;  void  ratio; 
and  Westergaard's  modulus  of  subgrade  reaction),  which  are  derived  from 
a variety  of  tests  (triaxial,  shear,  plate  bearing,  C3R,  and  penetrometer). 
Using  this  array  of  date,  materials  are  either  linearized  by  computation 
of  appropriate  constants,  or  a nonlinear  time-dependent  mathematical 
model  is  constructed  from  the  material  data  and  a behavioral  hypothesis. 

Any  practical  application  of  nonlinear  time-dependent  behavior  must  be 
done  withir.  a numerical  procedure.  Further  discussion  of  constitutive 
laws  is  presented  in  Section  3. 

2.1.fc  BOUNDARY  CONDITIONS 

Prescribing  the  boundary  conditions  for  an  idealization  is 
generally  broken  into  two  parts:  displacement  and  force.  For  the 

"principal"  techniques  this  requires  prescription  of  forces  which 
simulate  aircraft  and  gravitational  loads.  Additionally,  because  of 
the  "soil  island"  nature  of  finite  element  schemes,  displacements  that 
attempt  to  simulate  the  actual  continuous  nature  of  soil  must  be  pre- 
scribed along  the  edges  of  the  idealization. 

2.2  MATHEMATICS  ASSOCIATED  WITH  THE 
VESTERGAAFD  IDEALIZATION 

The  Weste>-gaard  idealization  transforms  the  airfield  problem  to 

the  oroblem  of  a surface  load  applied  to  a plate  on  a Winkler  foundation 

2 

(Figure  1).  As  normally  used,  the  method  considers  the  subgrade  to 
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have  a stiffness  k ; the  pavement  to  be  characterized  by  Young's 
modulus  E , Poisson's  ratio  v , and  pavement  thickness  h ; and  the 
load  to  be  circular  and  of  uniform  pressure  P and  radius  a . For 
this  situation  the  peak  deflection  is: 


6 = 


(1) 


where  £ is  the  radius  of  relative  stiffness  and 


^ = Eh3 

12(1  - v2)k 


and  the  maximum  tensile  stress  in  the  pavement  is: 


a 


3P 

8irh2 


(l  + v)  log 


(2) 


These  equations  are  a first-term  approximation  to  the  Westergaard 

idealization , and  are  generally  adequate  for  computing  the  idealization's 

response  when  the  top  layer  is  very  stiff  compared  to  the  other  layers. 

The  Westergaard  first-term  approximate  solutions  are  also  available  for 

idealizations  of  loads  applied  at  pavement  edges  and  joints,  and  for 

2 

elliptic  tire  prints. 

More  rigorous  solutions  for  the  Westergaard  idealization  may  be 
found  in  References  3 and  U.  Reference  U also  presents  solutions  for 
this  idealization  based  on  the  mathematics  for  the  elastic  layer 
idealization.  Since  the  Westergaard  idealization  is  a rather  crude 
approximation  of  the  real  problem,  more  rigorous  solutions  of  this 
idealization  are  seldom  warranted.  Comparisons  of  Westergaard  and 
elastic  layer  analyses  are  presented  in  Section  U.l. 


2.3  MATHEMATICS  ASSOCIATED  WITH  THE 
ELASTIC  LAYER  IDEALIZATION 


The  elastic  layer  idealization  transforms  the  airfield  problem 
to  the  problem  of  a load  applied  to  an  elastic,  horizontally  uniform, 
layered  system  (Figure  2).  As  normally  used,  the  method  considers  N 


bonded  layers  where  each  layer  is  defined  by  E , v , and  h . Each 
layer  extends  horizontally  to  infinity  with  the  bottom  layer  also  ex- 
tending vertically  to  infinity.  The  applied  load  of  uniform  pressure  P 
is  circular  with  radius  a . This  technique  requires  solution  of  rather 
complex  partial  differential  equations  which  are  usually  solved  numeri- 
cally with  the  assistance  of  a computer.  To  compute  responses  similar 
to  the  Westergaard  equations  (Equations  1 and  2)  requires  solution  of 
the  following  integral  equations  for  peak  deflection:* 


Pa(l  + v ) f 

5 = J K 

o 

J,  (x) 

+ ■ 2^a3  ~ ~ — dx 


(3) 


and  for  maximum  tensile  stress  in  the  top  layer: 


a2  + '“6  + al+^1*vi  + 


-ph 

+ ph^Ug  - au)]  e lj^x)  dx 


(U) 


The  a constants  are  determined  ^rom  boundary  conditions  while  is 

a Bessel  function,  p is  a transform  parameter,  and  x is  the  product 
of  p and  a . The  derivation  for  these  expressions  is  given  in 
Reference  h . 

Several  conputer  programs  are  available  for  computation  of  elastic 
layer  responses  (References  b-7).  Most  programs  allow  consideration  of 
either  bonded  or  frictionless  layer  interfaces;  results  of  such  considera- 
tions are  reported  in  Reference  8.  All  programs  require  numerical 
quadrature  of  integral  equations  and  evaluation  of  Bessel  functions; 

2 Rad  - v2) 


* For  a one-layer  system  (the  Boussinesq  problem),  6 


some  programs  perform  these  functions  more  accurately  than  others.  The 
results  obtained  from  this  idealization  for  linear,  elastic,  horizontally 
uniform,  layered  pavement  systems  are  "exact"  within  the  context  of  the 
theory  of  elasticity,  assuming  numerical  errors  are  negligible.  Thus, 
elastic  layer  responses  may  be  used  to  check  other  linear  elastic  tech- 
niques. Section  U presents  comparisons  of  elastic  layer  responses  with 
those  of  Westergaard  and  finite  element,  as  well  as  comparisons  of  com- 
putations derived  from  various  elastic  layer  programs. 

2.U  MATHEMATICS  ASSOCIATED  WITH  THE 
FINITE  ELEMENT  IDEALIZATION 

In  the  finite  element  analysis  of  airfield  pavements,  the  con- 
tinuous pavement  system  is  idealized  by  a finite  number  of  elements 
that  are  related  to  one  another  through  common  points,  called  nodes 
(Figure  3).  Each  element  and  node  are  assigned  a number,  and  the  bulk 
of  the  data  preparation  consists  of  inputting  to  the  computer  the  coordi- 
nates of  each  node  point  and  the  node  numbers  and  material  type  associated 
with  each  element.  Additional  data  concerning  material  characterization 
(e.g.,  E and  v)  for  each  material  type,  locations  of  loads,  and 
displacement  boundary  conditions  are  also  required.  These  data  allow 
the  computer  to  form  a mathematical  idealization  of  the  pavement  system. 
The  accuracy  of  this  idealization  is  dependent  on  a wide  variety  of 
factors  including  the  type  of  spatial  approximation  and  material  charac- 
terization used,  the  experience  of  the  engineer  in  preparation  and 
interpretation  of  the  results*  and  the  type  of  computer  hardware  avail- 
able. An  effort  will  be  made  throughout  the  report  to  elucidate  these 
factors. 

While  it  is  possible  to  obtain  a reasonably  good  working  knowl- 
edge of  classical  techniques  through  a limited  number  of  references 
(e.g.,  2-5),  obtaining  a good  working  knowledge  of  finite  element 
techniques  requires  familiarity  with  significantly  more  references  and 
considerable  opera* *<■>". el  »vp*ri»»nc#  with  the  finite  element  programs 
themselves.  The  authors  have  found  References  9-12  useful  in  providing 
basic  knowledge  of  the  finite  element  process.  More  detailed  knowledge 
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Figure  3.  Pavement  idealization  using  the  finite  element  procedure 

with  respect  to  the  basic  procedure  and  airfield  analysis  is  available 
in  various  reports  and  articles  which  are  mentioned  throughout  this 


report.  Unfortunately,  all  the  material  taken  as  a whole  does  not  re- 
place experience  or  the  necessity  for  good  "structural"  sense.  Even  the 
simplest  finite  element  idealizations  are  relatively  easy  to  do  poorly 
with  respect  to  accuracy  of  results,  efficient  utilization  of  computer 
hardware,  and  magnitude  of  engineer /programmer  effort  required.  While 
this  situation  is  certainly  cause  for  concern,  it  is  possible  to  obtain 
finite  element  results  easily  that  either  closely  duplicate  those  of 
classical  techniques  or  that  provide,  to  a high  degree  of  confidence, 
responses  that  only  the  finite  element  procedure  can  compute. 

2.U.1  SELECTION  OF  A FINITE 
ELEMENT  IDEALIZATION 

i Four  classes  of  finite  elenent  idealization  are  aDtO  icable  to  the 

analysis  of  pavement  systems.*  These  classes  are  grouped  recording  to 

* Within  each  class  many  permutations  of  the  finite  element  procedure 
are  permissible,  e.g.,  consideration  of  linear  and/or  nonlinear  materials. 
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spatial  limitations  on  the  variables  Involved  and  are:  (a)  plane  strain 

slice,  (b)  axisymmetric  solid,  (c)  prismatic  solid,  and  (d)  three- 
dimensional  solid.  A brief  description  of  each  class  follows: 

a.  A plane  strain  idealization  considers  a slice  of  the  pave- 
ment in  which  the  structure  is  limited  to  two  dimensions 
(Figure  M.  Only  results  within  the  slice  are  computed, 

and  changes  in  load  or  material  along  the  Z-axis  are  neglected. 
Because  generally  more  appropriate  classes  of  idealizations, 
i.e.,  (b),  (c),  and  (d),  have  become  available,  usage  of 
plane  strain  analysis  outside  research  activities  is  almost 
nonexistent.  A joint  analysis  using  this  class  is  shown  in 
Section  U.3. 

b.  The  axisymmetric  solid  idealization  considers  a layered, 
solid,  conical  frustrum  or  cylinder  of  pavement  loaded 
circularly  at  its  pole  (Figure  5).  A variation  of  the 
method  is  achieved  by  superimposing  results  of  one  analysis 
upon  another  to  simulate  a multiwheel  load  (Figure  6).  This 
type  of  idealization  has  been  widely  employed  in  a number 

of  different  computer  codes  for  pavement  analysis  (e.g., 
References  1,  13-17).*  This  technique  usually  requires  the 
least  cost  and  effort  to  use,  but  it  lacks  the  capability  to 
idealize  other  than  layer  systems  and  cannot  perform  nonlinear 
analysis  of  multiwheel  aircraft. 

c_.  The  prismatic  solid  idealization  considers  a nearly  general 
three-dimensional  system.  No  variation  of  geometrical  con- 
figuration is  allowed  from  cross  section  to  cross  section 
(X-Y  plane  in  Figure  7),  which  produces  a prismatic  solid 
along  the  Z-axis;  however,  loading  is  arbitrary.  This  class 
offers  a degree  of  structural  idealization  sophisticated 
enough  for  a variety  of  real  world  situations.  The  ease  of 
use  is  comparable  to  the  axisymmetric  analysis,  while  the 
cost  of  calculation  is  an  order  of  magnitude  greater. 
References  1 and  19-22  present  results  derived  from  the  AFPAV 
computer  program  for  prismatic  analysis. 

<1.  The  three-dimensional  idealization  considers  a complete 
general  structure;  material  and  load  configurations  are 
arbitrary  (Figure  8).  The  difficulty  of  use  is  an  order  of 
magnitude  greater  than  the  other  methods,  and  the  cost  of 
calculation  is  two  orders  of  magnitude  greater  than  the 
axisymmetric  method.  Because  of  cost  limitations,  the  use 
of  this  idealization  is  limited  to  special  applications. 

One  of  these  is  described  in  Reference  23. 


• Most  of  these  codes  are  based  on  a general  axisymmetric  code  written 
by  E.  L.  Wilson,  Reference  18. 
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Figure  6.  Axisywnetric  idealization  of  multiwheel 
load  using  superposition 


Figure  7.  Prismatic  idealization 
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Figure  8.  Three-dimensional  idealization 


2. It. 2 MATHEMATICAL-  BASIS  OF  THE 
FINITE  ELEMENT  PROCEDURE 


The  following  brief  introduction  to  the  mathematics  associated 
with  finite  element  analysis  is  included  to  provide  context  for  sub- 
sequent discussions. 

The  theorem  of  minimum  potential  energy  is  used  to  derive  the 
governing  equations  for  the  finite  element  method  and  may  be  written 
as: 


5P  * 0 


(5) 


where 


■ 1 pi 


i*l 


(6) 


and 


>i  • | J{t)T  [C]  {e>  dV  - J(6}T  {t>  dS  - J 6}T  [f} 


dV 


(7) 


where 


I * number  of  elements 
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potential  energy  for  i element 
volume  of  element 

transformed  strain  tensor 

stress-strain  transformation  (i.e.,  the  constitutive 
equations ) 


,th 


,th 


strain  tensor 

surface  area  of  i”"  element 
transformed  displacement  field  for  the  i“"  element 
prescribed  surrace  tractions 
prescribed  body  forces  per  unit  volume 


To  evaluate  these  integrals  requires  that  within  an  element  the  possible 
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modes  of  deformation  be  restricted.  The  theorem  further  requires  that 
the  displacement  field  across  element  boundaries  be  continuous  and  that 
certain  basic  deformation  modes  be  present  in  all  elements  regardless 
of  class  of  idealization  or  element  type. 

An  acceptable  set  of  deformation  modes  for  a two-dimensional 
element  is  shown  in  Figure  9*  The  displacement  field  produced  from  the 
combination  of  these  modes  is: 


u = a.  * a_s  ♦ a t ♦ a.  st 
1 2 3 « 


V * 0^  ♦ OgS  ♦ a^t  + OgSt 


The  a terms  are  constants  which  determine  the  amount  of  participation 
for  each  mode,  s and  t are  the  element's  natural  coordinates,  and 
u and  v are  the  displacements  in  the  s and  t directions,  respec- 
tively. While  Equations  8 and  9 are  a convenient  way  to  determine  the 
types  of  deformation  to  which  the  element  is  limited,  they  are  not  con- 
venient expressions  for  incorporation  into  the  finite  element  solution 
procedure  in  that  the  generalized  constant  o has  no  physical  meaning. 

A more  appropriate  description  of  the  element's  displacement  field  is 
obtained  by  manipulating  Equations  8 and  9 so  that  the  element's  dis- 
placement field  (u,v)  is  expressed  as  a function  of  the  element's  nodal 

displacements  (u  ,v  ) and  interpolation  functions  (h  ).  Thus, 
n n n 

u ■ h^u^  ♦ h2u2  ♦ h3u3  ♦ h^u^  (10) 


v - h^  ♦ h2v2  ♦ h3v3  ♦ h^ 


where 


h,  ■ l/l»(l  ♦ s)(l  ♦ t) 
h0  ■ 1/i* (1  - s)(l  ♦ t) 
h3  - 1/U(1  - s ) (l  - t) 
h^  - l/l»(l  ♦ s ) ( 1 - t) 
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The  element  shapes  that  are  used  with  the  four  classes  of  finite 
element  idealizations  are  shown  in  Figure  10;  most  programs  use  elements1 
that  have  quadrilateral  cross  sections.  While  there  are  differences  in 
the  derivation  of  each  element  type  based  on  the  class  of  idealization, 
the  basic  deformation  modes  associated  with  an  element  type  are  the  same 
regardless  of  class.  The  selection  of  these  modes  and  the  derivation  of 
equations  similar  to  10  and  11  are  based  either  on  isoparametric  con- 
cepts or  on  an  assemblage  of  constant  or  linear  strain  triangles. 

Delving  into  element  technology  is  beyond  the  scope  of  this  report,  but 
it  is  important  to  know  that  there  are  many  ways  to  derive  the  basic 
element  stiffness  with  concomitantly  varying  results.  A number  of  basic 
element  types  are  shown  in  Figure  11.  References  10  and  2k  provide  a 
comprehensive  presentation  concerning  selection  of  deformation  modes. 

Using  the  classical  strain-displacement  relations,  the  strain 
tensor  is  defined  in  terms  of  nodal  displacement: 


{ c } = (SHU.) 


(12) 


where 


[B]  « strain  displacement  transformation 

{U.}  ■ nodal  displacements  tor  the  ith  element  (i.e.,  u , v ) 
l n n 

Combining  Fquations  7 and  12  and  taking  the  variation  results  in 
formation  of  the  element  stiffness: 

[K.Hu.J  - (Fi)  (13) 


where 


(Kj)  ■ the  element  stiffness 

(F^)  » forces  applied  tc  the  nodes  of  the  ith  element 


or 
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QUADRILATERAL  PLANE 


a.  PLANE  STRAIN 


b.  AXISYMMETRIC 


c.  PRISMATIC 


d.  THREE-DIMENSIONAL 


Figure  10.  Element  shapes 


a.  FOUR- NODE 
ISOPARAMETRIC 


b.  EIGHT- NODE 
ISOPARAMETRIC 


c.  FOUR-NODE  QUADRILATERAL 
COMPOSED  OF  FOUR  CST 


d.  FOUR -NODE  QUADRILATERAL 
COMPOSED  OF  TWO  LINEARLY 
CONSTRAINED  LST 


LEGEND 

CST  CONSTANT  STRAIN  TRIANGLE 
LST  LINEAR  STRAIN  TRIANGLE 
INTERNAL  ELEMENT  EDGE 
O NOOE  POINT 

□ INTERNAL  DEGREE  OF  FREEDOM 


e.  EIGHT- NODE  QUADRILATERAL 
COMPOSED  OF  TWO  LST 


Figure  11.  Standard  element  typer. 
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[Kj  = J [B]T  [C]  [B]  dV 
V 

Adding  together  [K^  and  {Fi > for  all  I elements  results  in  the 
formation  of  the  relation  between  nodal  forces  and  displacements  for 
the  whole  pavement  structure.  That  is,  a set  of  algebraic  equations 
is  derived: 


[K] {U}  = {F}  (Ik) 

where  [K]  is  the  stiffness  for  the  whole  pavement  structure.  The 
unknowns,  {U}  , are  the  nodal  displacements  while  the  constants,  {F}  , 
are  the  prescribed  loads  applied  to  each  nodal  point.  After  displace- 
ment boundary  conditions  are  prescribed,  the  set  of  equations  is  solved 
using  standard  linear  algebraic  solution  techniques.  This  yields  values 
of  displacement  for  each  nodal  point  (i.e.,  {U}).  By  applying  Equa- 

tion 12  for  each  element,  the  element  strains  are  calculated.  The 
stresses  are  obtained  from  the  strains. 

For  a linear  structural  idealization,  [K]  remains  constant 
irrespective  of  the  magnitude  of  the  idealization's  displacement,  strain, 
or  stress.  The  nonlinear  analysis  of  a structure  primarily  implies 
that  [K]  changes  in  some  manner  as  a function  of  displacement,  stress, 
or  strain.  There  are  two  widespread  methods  employed  to  perform 
nonlinear  analysis.  One  is  known  as  the  initial  strain  method;  the  other 
is  the  tangent  method.  From  a user  viewpoint  these  two  methods  are  not 
very  different,  though  with  reference  to  the  computer  operations,  signif- 
icant differences  in  efficiency  and  approach  occur.  Both  methods  are 
appropriate  for  pavement  analysis  and  are  used  not  only  to  model  nonlinear 
materials,  but  to  utilize  large  deformation  theory,  to  simulate  material 
cracking,  and  to  allow  debonding  at  element  interfaces.  Discussion  of 
the  merits  of  these  methods  in  particular  situations  is  beyond  the  scope 
of  this  report.  References  25-27  are  suggested  for  further  study. 
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2.4.3  SELECTION  OF  A MESH 

Selection  of  a mesh  encompasses  a variety  of  decisions:,  choice 
of  element  sizings  and  spacings,  determination  of  optimal  pattern  for 
node  and  element  numbering,  selection  of  the  mesh  domain  (i.e.,  size  of 
soil  island),  and  prescription  of  the  boundary  conditions.'  The  mesh 
selected  depends  on  the  class  of  idealization  (e.g.,  prismatic  or 
axi symmetric)  and  the  permutation  involved  (e.g.,  nonlinear  or  linear 
axisymmetric ) ; the  type  of  element  employed  (e.g.,  constant  strain 
triangle  or  linear  strain  triangle);  and  the  type  of  structure  (e.g., 
flexible  versus  rigid  pavement).  Also  important  are  the  type  of  computer 
hardware  available  and  the  level  of  effort  to  be  expended.  Fortunately, 
for  most  applications,  only  the  axisymmetric  and  prismatic  solids  are 
of  practical  interest;  and  most  computer  codes  available  in  these  classes 
use  elements  that  have  essentially  linear  displacement  fields  (shown  in 
Figure  11a  and  11c).  Using  these  element  types  results  in  the  stress- 
strain  fields  that  are  nearly  constant  within  an  element.  This  result 
is  the  primary  factor  in  determining  mesh  selection:  where  the  gradient 

of  the  stress-strain  field  is  large,  more  elements  are  needed;  and  con- 
versely, where  the  stress-strain  is  relatively  constant,  fewer  elements 
are  needed.  The  mesh  is  also  affected  by  geometry  and  material 
boundaries. 

Figure  12  depicts  an  axisymmetric  finite  element  model  of  a re- 
paired bomb  crater  in  an  airfield  pavement;  this  study  is  reported  in 
Reference  17.  Small  elements  are  placed  directly  under  the  tire  where 
the  stress-strain  gradient  is  high.  Large  elements  are  used  farther  from 
the  load  because  the  stress-strain  field  is  relatively  constant.  Various 
odd-shaped  elements  are  used  to  match  the  material  boundaries.  The 
outside  edge  of  the  mesh  is  slanted  outward  to  accommodate  the  spreading 
of  the  load.  The  depth  of  the  mesh  is  fixed  at  a depth  where  no  signif- 
icant strain  occurs  in  the  soil.  For  this  case,  a trial-and-error 
scheme  was  employed  in  which  element  layers  were  successfully  added  to 
the  bottom  of  the  mesh  (thus  increasing  its  depth)  until  the  change  in 
peak  displacement  between  successive  runs  was  less  than  10  percent.  The 
width  of  the  mesh  is  set  by  a similar  process  while  the  slope  of  the  side 
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is  usually  taken  as  30  deg,*  based  on  results  reported  in  Reference  28. 
The  displacement  boundary  conditions  chosen  were  fixed  along  the  bottom 
and  free  or  horizontally  fixed  along  the  sides. 

An  additional  criterion  in  mesh  selection  is  ease  of  data  prepara- 
tion. Figure  12b  represents  an  alternative  mesh  to  that  of  Figure  12a. 

To  minimize  preparation  of  element  and  nodal  data,  a regular  lattice 
mesh  has  been  constructed.  The  material  interfaces  within  the  soil  mass 
are  only  approximately  defined  by  changing  the  element  material  numbers 
at  the  appropriate  places.  Justification  for  this  sacrifice  of  detail 
lies  in  the  approximate  nature  of  the  crater's  definition.  To  match 
this  approximation  exactly,  as  in  Figure  12a,  is  unwarranted . 

These  criteria  present  a reasonable  set  of  specifications  for  this 
pavement  model;  however,  under  unusual  circumstances  or  for  other  types 
of  models,  different  specifications  may  be  necessary. 

The  difficulties  most  often  associated  with  finite  element 
analysis  result  from  an  engineer's  "getting  in  over  his  head."  In- 
variably an  engineer  will  proceed  directly  to  solve  his  problem,  re- 
gardless of  its  nature.  Complicated  problems  must  be  solved  in  stages; 
experience  Is  extremely  valuable  in  designing  a competent  mesh  and 
producing  sav is factory  results. 

In  the  authors'  opinion,  the  best  procedure  for  acquiring  this 
experience  consists  of  starting  with  a simple  linear  problem  (e.g.,  that 
shown  in  Figure  13)  and  varying  each  mesh  parameter  (e.g.,  element 
distribution  and  number,  depth  and  width  of  mesh).  In  some  circumstances 
these  variations  will  produce  significant  changes  in  response.  The 
accuracy  of  results  is  assessed  by  comparison  with  an  elastic  layer 
solution  which  is  "exact"  for  a linear,  layered  pavement.  For  the 
linear  system,  stress-strain  should  be  in  good  agreement  for  the  two 
prediction  techniques  while  finite  element  displacements  will  generally 
be  a few  percent  less  than  those  of  elastic  layer. 

* A table  of  factors  for  converting  U.  S.  customary  units  of  measurement 
to  metric  (SI)  units  is  given  on  page  7. 
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Figure  13.  Finite  element  model  of  pavement  system 

For  nonlinear  and  nonuniform-layer  systems,  a similar  procedure 
to  this  one  should  be  adopted.  From  Known  results  or  past  computations, 
the  system  parameters  should  be  gradually  changed  to  achieve  the  desired 

result . 
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3 MATERIAL  MODELS 

The  most  important  component  of  any  prediction  method  is  the  mate- 
rial model  utilized  to  obtain  the  structural  responses.  It  is  para- 
doxical that  the  material  model  is  also  the  most  elusive  and  misused 
relationship  in  pavement  problems.  That  is,  other  components  of  the 
pavement  problem  such  as  geometry,  boundary  conditions,  and  loading  are 
relatively  easy  to  define  accurately  compared  to  the  material  model. 

Much  of  the  difficulty  related  to  material  models  is  due  to  the  plethora 
of  proposed  methods  and  models  reported  in  the  literature.  Constitutive 
forms  ranging  from  linear  to  highly  nonlinear  and  time  dependent  have 
been  suggested.  In  short,  there  is  no  unanimity  of  opinion  on  material 
models  for  pavements.  Although  divergent  opinions  provide  an  atmosphere 
for  active  research,  they  do  little  to  help  the  plight  of  practicing 
engineers  who  desperately  need  guidelines  on  material  models  for  pavement 
analysis.  In  the  following  sections,  attempts  will  be  made  to  provide 
general  guidelines  for  selecting  material  models;  however,  it  must  be 
understood  that  there  is  no  universal  answer  to  the  problem.  Each  pave- 
ment system  is  a unique  problem;  and  the  material  model  must  be  chosen 
commensurate  with  *he  knowledge  of  the  materials,  the  availability  of 
analytical  tools,  and  the  specified  performance  criteria. 

Implementation  of  material  models  cannot  be  divorced  from  the 
analytical  technique  employed  for  the  pavement  analysis.  To  illustrate 
this  assertion,  a suggested  finite  element  framework  suitable  for  accom- 
modating most  material  models  is  presented  along  with  various  solution 
strategies.  After  this  framework  is  established,  material  models  are 
r | examined  and  discussed. 

[ I 

[ < 3.1  SUGGESTED  ANALYTICAL  FRAMEWORK 

FOR  FINITE  ELEMENT  CODES 

To  begin  with,  a finite  element  code  must  be  modularized  to  accom- 
modate new  material  laws  easily  as  they  become  available  or  otherwise  be 
doomed  to  extinction.  Secondly,  the  code  should  be  structured  for  in- 
cremental loading  as  opposed  to  one-step  total  loading  since  the  one-step 
leading  is  inherently  a subset  of  incremental  loading  and  most  all 


material  models  can  be  cast  in  an  incremental  form.  In  addition, 
algorithms  and  storage  for  iterating  within  the  load  step  should  be 
provided,  i.e.,  the  ability  to  apply  a load  step  repeatedly  until  both 
equilibrium  and  the  material  model  are  in  agreement.  These  notions  are 
described  in  a schematic  flow  chart  (Figure  lU). 

This  general  algorithm  allows  considerations  of  several  popular 
solution  strategies  for  treating  nonlinear  materials  including  the 
secant  method,  the  tangent  method,  the  modified  tangent  method,  and  the 
chord  method.*  Figure  15  illustrates  the  concepts  associated  with  each 
of  these  methods.  The  secant  method  implies  that  the  total  load  is 
applied  in  one  step  and  the  process  is  iterated  to  find  a secant  modulus 
satisfying  both  equilibrium  and  the  material  law  (stress-strain  curve). 
The  tangent  method  implies  that  the  load  is  applied  in  a series  of 
steps.  At  the  end  of  each  step  the  tangent  of  the  material  law  is 
evaluated  at  the  accumulated  stress-strain  level  to  provide  the  modulus 
for  the  next  load  step.  Note  that  the  stress-strain  responses  cal- 
culated by  this  method  increasingly  diverge  from  the  material  law  under 
monotonic  loading.  The  modified  tangent  method  avoids  this  divergence 
by  iterating  within  the  load  step  to  determine  a modulus  which  is  an 
average  of  the  material  law  tangent  at  the  beginning  and  end  of  the  load 
step.  The  chord  method  is  the  secant  method  applied  in  a step-by-step 
fashion;  the  chord  method  inscribes  the  material  law. 

The  selection  of  one  method  over  another  is  largely  a question 
of  cost  versus  accuracy.  The  secant  method  is  acceptably  accurate  in 
cases  of  near  proportional  loading  (i.e.,  principal  stresses  remain  in 
a near  constant  ratio  throughout  the  load  path)  and  is  generally  the 
least  costly  procedure.  The  tangent  method,  which  is  comparable  in  cost, 
is  more  appropriate  for  nonproportional  loading  and  also  provides  a 
history  of  the  response  of  the  structure.  Both  the  modified  tangent 
and  chord  methods  provide  improved  accuracy  over  the  secant  n..d  tangent 
methods;  however,  their  cost  is  significantly  greater. 

• These  solution  strategies  are  an  integral  part  or  (but  deferent  from) 
the  nonlinear  methods  employed  by  the  finite  element  prog  "m  mentioned 
at  the  end  of  Section  2.1*. 2. 
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Figure  15.  Solution  strategies  f->r  nonlinear  materials 


As  a general  rule  of  thumb,  the  authors  contend  that  all  things 
being  equal,  the  higher  cost  penalty  should  not  deter  engineers  from 
using  more  accurate  algorithms  since  the  computer  cost  is  usually  minis- 
cule compared  to  the  overall  project. 

3.2  AN  EXAMPLE  OF  THE  EMPLOYMENT  OF 
NONLINEAR  MATERIAL  MODELS  USING 
THE  FINITE  ELEMENT  PROCEDURE 

Analysis  of  nonlinear  materials  generally  consists  rf  a sequence 
of  "linearlike"  analyses  (referred  to  as  steps)  which  are  related  to 
one  another.  At  each  step  a linear  response  of  the  structure  (i.e., 
displacement,  stress,  and  strain)  is  computed.  The  material  parameters 
used  during  a step  for  each  element  are  determined  from  the  response  of 
the  element  computed  in  the  previous  step.  The  outcome  of  such  an 
analysis  is  a sequence  of  linear  responses.  Taken  as  a whole,  they 
represent  a nonlinear  response  to  the  applied  load  (mathematically,  this 
is  known  as  a linear  piecewise  approximation  of  a nonlinear  function). 

The  tangent  method  and  the  model  illustrated  in  Figure  13  will 
be  used  to  demonstrate  the  technique.  In  this  example  a single  wheel 
is  being  supported  by  a concrete  surface  upon  a single  layer  of  soil. 

To  simplify  the  discussion,  only  the  response  of  a single  soil  element 
will  be  considered,  marked  3 in  the  figure.  The  soil’s  nonlinear 
material  properties  are  shown  in  Figure  lo;  for  simplicity,  only  hydro- 
static stress  (i.e.,  (op  + ♦ sg)/3)  versus  volumetric  strain  (i.e., 

e„  ■ -(c0  +€_■♦•  ej)  will  be  considered  Material  unloading  and  pre- 
loading  are  omitted.  Five  steps  will  be  used  to  approximate  the  nonlinear 
response. 

For  each  step,  the  response  of  the  soil  element  will  be  computed 

for  one-fifth  of  the  total  load  (i.e.,  30  psi).  At  the  beginning  of 

the  ntJl  step,  the  element's  bulk  modulus  K is  computed  based  on 

n 

the  strain  in  element  S at  the  end  of  the  previous  step  en  . Thus, 
for  *he  first  step,  where  = 0 , is  the  tanaent  at  the  beginning 

of  the  material  curve.  Because  of  the  linear  approximation,  an  error 
is  accumulated  at  the  end  of  each  step  such  that  the  material  curve  the 
element  "sees"  shadows  the  curve  derived  in  the  laboratory.  The  degree 
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to  which  the  shadow  approximates  the  measured  curve  is  for  this  example 
improved  merely  by  increasing  the  number  of  steps  used  for  nonlinear 
analysis.  The  responses  computed  during  each  step  are  cumulatively 
added  to  those  of  the  previous  steps,  resulting  at  the  end  of  the  last 
step  in  the  response  of  the  pavement  to  the  total  applied  load  (for  this 
example,  100  psi). 

The  number  of  steps  used  for  the  analysis  is  a function  of  the 
degree  of  nonlinearity  present.  Five  steps  have  been  found  to  give 
reasonable  results  while  "«ing  more  than  20  steps  provides  little  added 
benefit.  The  procedures  for  this  study  were  to  use  five  steps  for  all 
the  preliminary  calculations  needed  to  establish  the  various  mesh  param- 
eters and  to  use  20  steps  for  the  final  runs.  (This  reduced  the  error 
associated  with  the  linear  piecewise  approximation  of  the  measured 
material  curve  to  less  than  10  percent.)  For  this  study,  only  where 
concrete  overlayed  a "soft"  soil  resulting  in  extensive  concrete  cracking 
were  there  significant  differences  between  the  5-  and  20-step  analyses. 

3.3  CLASSIFICATION  OF  MATERIAL 
MODELS  FOR  PAVEMENTS 

A convenient  method  for  categorizing  material  lavs  is  with  the 

descriptors:  time-dependent  ' ersus  time- independent  and  nonlinear  versus 

linear.  Here,  time-dependent  implies  that  a real-time  variable  is 

included  in  the  constitutive  relation,  and  nonlinear  implies  that  stress 

is  not  a linear  function  of  strain  when  all  other  state  variables  are 

held  constant.  These  descriptors  are  utilized  in  Figure  17  to  denote  a 

hierarchy  of  constitutive  foi  The  top  box  of  this  figure  represents 

the  most  general  form  (time-uependent  and  nonlinear)  of  which  examples 

are  viscoplasticity  and  nonlinear  creep  laws.  Much  work  has  been  done 

29 

on  these  constitutive  forms  for  metals  but  very  little  for  pevements. 

It  is  well  known  that  pavements  exhibit  significant  time-dependent 
behavior  particularly  for  saturated  coheuive  soils  and/or  asphaltic 
surfaces.  However,  unless  the  real-time  response  history  i*  of  interest, 
the  added  complication  of  using  time-dependent  models  is  questionable. 
That  is,  from  an  engineering  viewpoint  the  capacity  of  a pavement  system 


can  be  adequately  determined  by  considering  all  loads  applied  for  long 
durations,  thereby  accounting  for  creep  and  relaxation  effects.  No  doubt, 
future  research  may  make  these  time-dependent  nonlinear  forms  very  attrac- 
tive, particularly  if  failure  criteria  can  be  built  into  the  models. 
However,  until  such  time,  these  forms  belong  more  to  the  realm  of  research 
than  to  today's  design  tools. 

In  the  next  level  of  Figure  IT  the  general  constitutive  form  is 
subdivided  into  two  categories:  time-dependent  and  linear  and  time- 

independent  and  nonlinear.  Examples  of  the  first  category  include  visco- 
plasticity and  linear  creep  laws.  Again,  as  previously  discussed,  the 
current  benefits  from  including  real  time  in  the  constitutive  form  are 
marginal.  Worse  still,  the  implied  linear  stress-strain  assumption  is 
simply  not  in  accordance  with  observed  behavior  of  any  component  of  a 
pavement  system.  That  is,  all  pavement  components  that  exhibit  a sig- 
nificant time-dependent  response  also  exhibit  a significant  nonlinearity, 
e.g.,  saturated  granular  soils,  cohesive  soils,  and  asphaltic  cements. 
Consequently,  it  does  not  appear  that  time-dependent  and  linear  material 
models  will  ever  play  a significant  role  in  pavement  analysis  with  the 


Figure  17.  Classification  of  material  models 


possible  exception  of  consolidation  problems  or  pa -ements  composed  of 
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sea-ice. 

The  adjoining  category,  time-independent  and  nonlinear,  is  by  far 
and  away  the  most  significant  category;  consequently  this  category  will 
be  discussed  in  greater  depth. 

It  is  convenient  to  further  subdivide  the  time-independent  and 

nonlinear  classification  into  two  groups:  plasticity  models  and  variable 

modulus  models.  The  former  grouping  is  based  on  the  theory  of  plasticity 

which  in  general  requires  a yield  criterion,  a hardening  rule,  and  a 

flow  rule.  A yield  criterion  defines  the  onset  of  plastic  yielding  and 

is  usually  assumed  to  be  a function  of  the  stress  invariants.  The 

hardening  rule  redefines  the  yield  criterion  after  plastic  deformation 

has  occurred  and  is  usually  assumed  to  be  a function  of  plastic  work  and 

stress  level.  Lastly,  the  flow  rule  relates  increments  of  plastic  strain 

to  increments  of  stress  after  the  yield  criterion  is  satisfied.  Examples 

of  plastic  models  applied  to  pavements  are  the  Drucker-Prager , 
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Mohr-Coulomb , and  capped  models. 

From  an  academic  viewpoint,  the  plasticity  models  are  more 
glamorous  than  the  variable  modulus  models  (discussed  next)  because  they 
generally  satisfy  rigorous  theoretical  requirements  and  are  inherently 
capable  of  treating  unloading  and  cyclic  loading  for  fatigue  considera- 
tions. On  the  negative  side,  plasticity  models  generally  do  not  correlate 
well  with  triaxial  data  from  soil  specimens.  Also,  the  parameters  of 
the  plasticity  models  are  relatively  difficult  to  determine. 

Variable  modulus  models  are  based  on  the  hypothesis  that  stress 
increments  can  be  related  to  strain  increments  by  an  elastic  constitutive 
matrix  wherein  the  components  of  the  constitutive  matrix  are  dependent 
on  the  level  of  stress  and/or  strain  (i.e.,  (Ao)  - [C]{Ae}  where  (Ac) 
and  (Ae)  are  increments  of  stress  and  strain  and  [ C ] is  the  elastic 
constitutive  matrix  whose  components  are  dependent  on  the  current  total 
level  of  stress  and  strain).  The  variable  modulus  models  represent  mate- 
rials of  the  hypoelastic  classification;  that  is,  the  constitutive  com- 
ponents are  dependent  upon  initial  conditions  and  the  stress  path. 
Consequently,  the  term  nonlinear  elastic  is  not  appropriate  for  variable 


modulus  models  since  nonlinear  elastic  implies  path  independence.  The 
advantages  of  variable  modulus  models  are  their  inherent  ability  to 
approximate  experimental  data  closely  and  the  relative  ease  of  deter- 
mining the  parameters  of  the  model.  The  disadvantages  of  this  approach 
are  that  unloading  is  not  inherently  built  into  the  model  and  must  be 
treated  in  an  ad  hoc  fashion.  Also  special  programming  features  may 
have  to  be  employed  to  avoid  possible  numerical  and  theoretical  diffi- 
culties arising  from  energy  violations  (e.g.,  Poisson's  ratio  exceeding 
0.5).  Examples  of  variable  modulus  models  are  given  in  References  32-31*. 
Further  discussion  of  plasticity  and  variable  modulus  models  as  they 
relate  to  pavements  is  given  in  the  next  section. 

The  last  category  of  material  models  shown  in  Figure  17  is  the 
linear  form.  The  linear  form  is  a straightforward  application  of  the 
generalized  Hook's  law.  The  obvious  advantage  of  the  linear  form  is  its 
utter  simplicity;  it  also  permits  the  usage  of  classical  solution  tech- 
niques. Of  course,  the  disadvantage  of  the  linear  form  is  the  inability 
to  model  the  nonlinear  behavior  of  pavement  systems,  particularly  the 
base  course  and  subgrade  components. 

To  summarize,  it  is  felt  that  the  time-independent  and  nonlinear 
category  of  Figure  17  embodies  the  most  relevant  models  for  pavement 
analysis.  Time-dependent  and  nonlinear  models  generally  impose  a level 
of  sophistication  in  excess  of  today's  performance  criteria;  however, 
their  future  is  very  promising.  Time-dependent  and  linear  models  do 
not  now  or  in  the  future  appear  to  be  relevant  in  pavement  analysis. 
Lastly,  the  simple  linear  forms  have  limited  applications  and  may  be 
entirely  phased  out  in  future  developments. 


In  the  next  section,  particular  models  in  the  time-independent  and 
nonlinear  category  are  discussed. 

3.1*  MATERIAL  MODELS  FOR  PAVEMEUTS 

Wearing  courses  are  the  most  significant  component  of  a layered 
pavement  system  since  they  directly  receive  the  load  and  distribute  it 
to  the  layers  below;  proper  material  modeling  is  essential  for  meaningful 


results.*  Both  variable  modulus  and  plasticity  models  have  been  exten- 
sively used  for  these  materials  with  apparent  success.  In  particular , 
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the  hybrid  plasticity  model  developed  by  Kupfer  et  al.  is  well  suited 

for  concrete.  The  yield  surface  of  this  model  faithfully  matches  the 

data  points  of  concrete  specimens  in  all  quadrants  of  the  biaxial  stress 

plane.  Accordingly,  the  model  can  predict  the  onset  of  concrete  cracking 

in  tension  zones  as  well  as  nonlinear  responses  in  the  compression  zones. 

A finite  element  development  of  this  model  is  given  in  Reference  36. 

Base  course  layers  are  almost  exclusively  composed  of  various 

gradations  of  granular  materials  or  stabilized  soil  while  the  subgrade 

layers  may  include  granular,  mixed,  and  cohesive  soils.  A model  capable 

of  treating  any  of  these  soils  is  a variable  modulus  model  originated 
32 

by  Hardin.  The  objective  of  the  following  discussion  concerning  the 
Hardin  model  is  to  demonstrate  procedures  that  should  be  used  to  evaluate 
any  soil  model  and  to  outline  desirable  model  features. 

The  Hardin  model  presents  a hyperbolic  relationship  for  the  secant 
shear  modulus  as  follows: 

G 

„ _ max 

Gs-  — 

where 

G = secant  shear  modulus 
s 

G * maximum  shear  modulus,  a function  of  spherical  stress 
max 

Yh  * hyperbolic  shear  strain,  a function  of  shear  strain  and 
spherical  stress  for  a particular  soil 

The  beauty  of  Hardin's  approach  is  that  not  only  does  he  present 

the  form  of  the  relationships  for  G and  y.  , but  he  also  identifies 

max  n 

the  soil-dependent  parameters  and  provides  equations  for  evaluating  them 
in  terms  of  soil  type  (granular,  mixed,  or  cohesive),  void  ratio,  percent 

* An  error  in  material  characterization  for  the  top  layer  has  a dramatic 
impact  on  the  predicted  response  for  the  whole  pavement.  This  is  es- 
pecially true  for  asphalt  where  the  authors  have  observed  others  using 
values  of  E from  20,000  to  a million  psi,  and  for  concrete  where 
tensile  cracking  is  extensive. 
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Figure  l8.  G versus  y 
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* J.  T.  Ballard,  U.  S.  Army  Engineer  Waterways  Experiment  Station,  per- 
sonal communication  with  *he  Naval  Civil  Engineering  Laboratory, 

6 September  1973. 
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prediction  was  not  obtained  by  curve  fitting,  but  rather  by  a straight- 
forward application  of  the  model  where  the  only  soil  information  required 
was  void  ratio. 

This  result  implies  that  it  is  possible  to  utilize  the  Hardin 
shear  model  in  the  absence  of  test  data,  a situation  often  encountered 
by  pavement  analysts. 

By  itself,  Hardin’s  relationship  for  the  shear  modulus  does  not 
provide  a complete  variable  modulus  model.  To  complete  the  constitutive 
form,  it  is  necessary  to  incorporate  at  least  one  additional  elastic 
parameter  (function)  into  the  constitutive  matrix.  For  example,  Young's 
modulus,  bulk  modulus,  or  Poisson's  ratio  could  be  selected.  Of  these 
choices,  Poisson's  ratio  is  the  most  convenient  based  on  the  consideration 
of  avoiding  numerical  difficulties  and  maintaining  theoretical  energy 
requirements.  By  developing  a function  for  Poisson's  ratio  such  that  its 
range  is  within  the  limits  0.0  1 v < 0.5  , the  energy  requirements  are 
satisfied  regardless  of  the  current  value  of  the  shear  modulus.  On  the 
other  hand,  if  the  bulk  modulus  or  Young's  modulus  were  used,  their  valid 
range  is  dependent  on  the  current  value  of  the  shear  modulus  (i.e.,  bulk 
modulus  must  be  greater  than  two-thirds  of  the  shear  modulus,  and  Young's 
modulus  must  be  greater  than  two  but  less  than  three  times  the  value  of 
the  shear  modulus). 

Many  investigators  assume  Poisson's  ratio  is  constant  in  soil 
models;  however,  this  assumption  is  based  more  on  convenience  than  on 
actual  material  behavior.  The  fact  that  Poisson's  ratio  does  vary 
during  loading  is  demonstrated  in  Figure  19-  Here  the  observed  Poisson's 
ratio  of  Cook's  bayou  sand*  is  plotted  for  four  triaxial  tests  as  a 
function  of  normalized  shear  strain.  Clearly  it  is  observed  that 
Poisson's  ratio  is  a function  of  both  shear  strain  and  confining  pressure 
and  varies  from  approximately  0.1  to  0.5.  Accordingly,  any  serious  at- 
tempt to  model  soil  behavior  must  incorporate  a relationship  for  Poisson's 
ratio.  To  this  end,  a general  relationship  for  Poisson's  ratio  has  been 
developed  in  a manner  analogous  to  Hardin's  shear  modulus.  The  interested 
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Figure  19.  Poisson's  ratio  v versus  ratio  of  measured  shear 
strain  y to  a reference  shear  strain  Yr 

reader  may  find  a complete  presentation  of  this  development  in  Refer- 
ence 37  along  with  a review  of  Hardin's  and  other  relationships. 

3.5  SUMMARY 

The  main  points  of  this  section  are  summarised  below: 

a.  Finite  element  codes  should  be  structured  to  operate  on  an 
~ incremental  basis  with  allowance  for  iterating  within  the 

load  step. 

b.  Material  models  should  be  isolated  in  subroutines  to  permit 
~ easy  additions,  deletions,  and/or  modifications. 

c.  Solution  strategies  for  nonlinear  equations  should  not  be 

~ compromised  to  promote  efficiency  at  the  expense  of  accuracy. 

d.  Time-dependent  and  nonlinear  material  models  are  for  the  most 
part  overly  sophisticated  with  respect  to  the  currently  em- 
ployed performance  criteria. 

e.  Time-dependent  and  linear  material  models  have  little  relevance 
in  pavement  analysis. 
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f . Time- independent  and  nonlinear  material  models  are  most  rele- 
vant to  the  pavement  problem.  In  particular,  plasticity 
models  are  adequate  for  the  wearing  courses,  and  variable 
modulus  models  are  adequate  for  the  soil  layers. 

Hardin's  shear  modulus  relationship  is  an  excellent  example 
of  the  generality  and  versatility  that  should  be  incorporated 
into  soil  models. 


4 COMPARISON  AND  PRESENTATION  OF  ANALYTICAL 
AND  MEASURED  RESULTS 


Previous  sections  have  dealt  in  a general  way  with  particular 
facets  of  prediction  techniques.  This  section  presents  severed  applica- 
tions of  prediction  techniques  and  by  illustrating  their  application  and 
results  ties  them  to  the  foregoing  general  discussion.  The  techniques 
covered  include  Westergaard,  elastic  layer  (Burmister  problem),  and 
linear  and  noalinear  finite  element.  Most  of  the  examples  presented 
have  been  taken  from  various  previously  published  papers.  Comparison 
of  one  technique  versus  another  is  shown  wherever  available. 

4.1  WESTERGAARD  AND  ELASTIC  LAYER 
ANALYSES 

Figure  20  shows  various  pavement  sections  which  are  used  to 
demonstrate  the  application  of  Westergaard  and  elastic  layer  analysis 
techniques.  The  first  five  are  flexible  sections  while  the  second  five 
are  rigid.  Elastic  layer  responses  are  shown  in  Table  1 for  the  sec- 
tions of  Figure  20  using  four  different  elastic  layer  computer  codes: 
ELAST  (4),  BISTRO  (5),  CRANLAY  (7),  and  CHEVRON  (6).  Because  these 
sections  are  assumed  to  be  linear  elastic  and  have  no  Joints  or  edges, 
the  layer  elastic  response  is  "exact"  within  the  context  of  the  theory 
of  elasticity,  assuming  negligible  computational  errors. 

The  Westergaard  equations  shown  in  Section  2.2,  which  are  those 
standardly  used,  are  a first  approximation  for  the  solution  of  a plate 
on  a fluid  foundation  (Figure  1).  Westergaard  in  Reference  3 indicates 
that  including  more  than  the  first  term  of  the  Kei  and  Ker  function 
yields  more  accurate  results.  Responses  for  the  single  term  (Equa- 
tions 1 and  2)  and  a six-term  approximation  are  shown  in  Table  1.  Both 
of  these  approximations  assume  that  the  top  layer  contains  a neutral 
axis;  erroneous  results  occur  when  this  is  untrue.  To  eliminate  this 
problem,  the  Westergaard  idealization  is  solved  using  mathematics  similar 

to  that  used  for  the  layer  elastic  idealization.*4  in  this  case,  the 
first  layer  is  considered  an  axisymmetric  solid  rather  than  a plate. 

The  seventh  column  of  Table  1 contains  these  results. 
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Figure  20.  Selected  pavement  sections 


The  Westergaard  and  the  elastic  layer  idealizations  may  be  re- 
lated through  their  common  mechanistic  basis  using  the  following  function 
(see  Figures  1 and  2 for  symbols). 

k = F(h2E2v2,  h3E3v3...ENvN) 

The  derivation  of  F is  presented  in  detail  in  Reference  U,  and  the 
results  of  its  usage  are  shown  in  Table  1 as  kw  . 

To  enable  direct  comparison  of  elastic  layer  and  Westergaard 
results,  all  of  the  Westergaard  responses  shown  in  the  table  are  com- 
puted using  . Thus,  at  least  for  the  sections  of  Figure  20,  the 
predictions  of  the  Westergaard  idealization  (whichever  solution  technique 
is  used,  i.e.,  columns  5,  6,  or  7 of  Table  l)  are  "spotty"  when  measured 
against  the  "exact"  elastic  layer  solution  (computed  by  any  of  the  four 
codes  shown). 

Values  for  the  Westergaard  subgrade  modulus  k are  shown  where 
available  in  the  last  column  of  Table  1.  Comparison  of  ktf  (computed 
from  layer  properties)  with  k produces  a measure  of  the  accuracy  with 
which  the  elastic  layer  material  properties  predict  the  results  of  a 
plate  bearing  test.  In  this  manner  k^  can  be  used  to  assess  the  over- 
all performance  of  the  elastic  layer  parameters.  For  the  sections  of 
Figure  20  where  most  of  the  E and  v values  were  based  on  rather 
sketchy  data,  it  is  expected  that  k and  ky  will  not  agree  very  often. 

The  advantages  of  both  these  techniques,  especially  Westergaard, 
are  their  ease  of  application  and  relatively  modest  requirement  for 
preparation  of  data  or  interpretation  of  results.  If  these  techniques 
are  carefully  applied  with  knowledge  of  their  limitations,  useful  design 
data  can  be  generated.  As  an  analysis  procedure,  however,  these  tech- 
niques are  at  a distinct  disadvantage  with  respect  to  all  but  linear 
layered  pavements.  The  crux  of  the  matter  is  that  while,  for  example, 
a nonlinear  k can  be  used  for  a Westergaard  analysis,  the  lack  of 
capacity  in  these  formulations  to  assess  what  these  nonlinear  approxima- 
tions mean  makes  such  results  suspect.  Possibly  use  of  more  refined 
methods,  such  as  finite  element,  to  prove  the  validity  of  such  extreme 
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extensions  of  relatively  simple  theories  vill  make  such  approaches 
desirable. 

4.2  LINEAR  FINITE  ELEMENT 

Finite  element  techniques  have  two  significant  differences  from 
the  linear  idealizations  of  Section  4.1:  the  capacity  to  vary  material 

properties  within  the  layer  and  the  consideration  of  a soil  island. 

A single-  and  multiwheel  analysis  for  Section  4 of  Figure  20  is  presented 
to  establish  the  relationship  between  finite  element  and  elastic  layer 
analyses  and  to  study  the  peculiarities  of  the  finite  element  procedure. 
This  pavement  section  is  Item  4 of  the  C-5A  test  track  at  Waterways 
Experiment  Station  (WES). 

4.2.1  SINGLE  WHEEL  ON 
LAYERED  PAVEMENT 

Figure  21*  shows  the  response  predicted  by  an  axisynsnetric  finite 
element  code  (WIL67)  and  an  elastic  layer  code  (BISTRO)  for  a single 
C-5A  wheel  loaded  by  30,000  lb.  Several  features  are  illustrated  by 
this  figure: 

a.  The  stress  responses  predicted  by  elastic  layer  and  finite 

— element  codes  are  identical  while  the  deflections  are  dif- 
ferent by  a constant  amount.  This  is  a typical  result  where 
the  disparity  is  attributable  to  the  soil  island  nature  of 
the  finite  elemert  procedure;  that  is,  the  values  of  vertical 
deflection  are  dependent  on  the  integration  of  all  the  vertical 
strains  along  a vertical  line.  If  the  mesh  is  truncated 
before  these  strains  are  negligible,  then  a portion  of  dis- 
placement is  lost,  which  leads  to  a constant  error  factor 
between  the  "exact"  elastic  layer  and  finite  element  idealiza- 
tions. However,  the  stresses,  strains,  and  deflection 
shapes,  which  are  derivatives  of  the  vertical  displacement, 
are  easily  matched  between  the  two  methods  for  any  reasonably 
designed  mesh.  To  match  the  deflection  requires  a convergence 
study  which  is  discussed  in  Section  4.2.3. 

b.  The  linear  pavement  response  is  relatively  insensitive  to  the 

~ distribution  of  stiffness  in  the  lower  layers.  This  is  demon- 
strated by  comparison  of  the  responses  for  the  seven  actual 
layers  to  a smeared  four-layer  system.  The  smeared  values  are 
shown  in  Figure  21  as  E£  and  v£  . 

• This  figure  and  the  cost  data  were  extracted  from  Reference  1. 
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c.  The  reversed  curvature  form  of  the  measured  displacements  is 
typical  of  test  results  and  is  indicative  of  nonlinear 
behavior.  The  rather  uniform  curvature  of  the  predicted 
displacements  is  typical  of  linear  layered  analysis,  irrespec- 
tive of  the  method  applied. 

d.  Stiffness  within  a layer  must  be  nonuniform  to  match  the 

test  results.  The  disparity  between  measured  and  predicted 
stress  (Figure  21a)  can  be  attributed  to  two  possible 
failings  of  the  model:  first,  that  the  material  of  layer  4 

is  too  flexible,  or  second,  that  too  much  flexural  stiff- 
ness is  present  in  layer  3.  Examination  of  these  two  pos- 
sibilities yields  the  following:  if  the  flexural  stiffness 

E3  is  decreased,  the  magnitude  and  shape  of  the  vertical 
deflection  curve  will  become  larger  and  less  vertical  in 
the  third  layer  than  that  shown  in  Figure  21b.  Second,  if 
the  fourth  layer's  stiffness  E^  is  increased,  the  deflected 
shape  across  the  fourth  layer  becomes  steeper.  Neither  of 
these  courses  will  produce  an  improvement  in  the  predicted 
responses,  and  moreover,  little  change  will  occur  in  the 
stress  field  because  of  equilibrium  considerations  associated 
with  layered  systems.  Thus,  neither  possibility  offers  the 
answer,  and  in  fact,  this  disparity  illustrates  the  funda- 
mental flaw  associated  with  elastic  layer  idealizations: 
that  is,  uniformity  of  stiffness  within  a layer.  The  only 
reasonable  way  that  the  measured  responses  for  this  section 
can  be  matched  by  analytical  prediction  is  by  resorting  to 

a nonuniform  layer  stiffness,  that  is,  (l)  by  an  increase 
of  stiffness  in  layer  3 directly  under  the  load  which  causes 
an  increase  in  vertical  stress,  but  little  increase  in 
flexural  stiffness;  and  (2)  by  a decrease  in  stiffness  at 
the  center  of  layers  4 and  5 under  the  load  to  provide  the 
reverse  curvature  of  the  WES  test  results. 

The  computer  time  for  running  the  BISTRO  code  is  proportional  to  the 
number  of  locations  at  which  responses  are  requested  and  the  number  of 
pavement  layers  (for  Figure  21,  4 sec  per  location  on  a CDC  6600  for 
four  layers  and  6 sec  for  seven  layers).  For  the  WIL67  code  the  time 
is  related  to  the  mesh,  which  is  selected  relatively  independent  of  the 
number  of  responses  desired  and  pavement  layers  (4o  sec  on  a CDC  6600 
for  a 480-element  mesh).  Thus,  if  more  than  10  responses  are  to  be  com- 
puted, use  of  the  finite  element  technique  would  be  cheaper  in  this  case. 
Because  of  the  data  generators  associated  with  WIL67,  the  amount  of  data 
preparation  is  the  same  for  both  codes. 
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4.2.2  MULTIPLE-WHEEL  GEAR  ON 
LAYERED  PAVEMENT 

While  single-vheel  analyses  have  much  edacational  value  and  are 
of  practical  significance  at  installations  predominately  used  by  fighter 
aircraft,  most  airfields  are  designed  for  the  large  multiwheeled  aircraft 
such  as  the  C-5A  shown  in  Figure  22.  To  obtain  multiwheel  results  re- 
quires use  of  superposition  for  elastic  layer  and  axisymmetric  idealiza- 
tions, or  results  may  be  directly  obtained  using  prismatic  or  three- 
dimensional  solid  idealizations.  Figure  23  shows  the  predicted  response 
for  the  12-wheel  assembly  placed  on  the  four-layer  approximation  of 
Section  4.  Results  similar  to  those  of  Figure  21  are  obtained. 

&.  Figure  23a  shows  the  predicted  deflection  basin  along  the  gear 
centroid  (i.e.,  along  the  X-axis).  Two  different  mesh  depths 
were  taken,  both  fairly  shallow,  causing  considerable  disparity 
between  finite  element  and  elastic  layer. 

b.  The  data  presented  in  Figure  23b  and  c show  the  maximum  re- 
sponse recorded  within  a horizontal  plane  as  a function  of 
depth;  that  is,  these  curves  are  envelopes  of  peak  response 
(denoted  limiting  curves  in  Figure  23).  Also  shown  in  these 
curves  is  the  effect  of  increasing  the  stiffness  of  layer  4: 
no  effect  occurs  in  the  stress;  and  while  the  displaced 
shape  is  the  same,  the  curve  is  shifted  to  the  left. 

£.  Figure  23d  and  e shows  the  horizontal  and  shear  stress  under 
the  points  shown  in  Figure  22.  These  values  of  stress  are 
also  unaffected  by  changing  E^  to  8000  psi.1 

The  time  for  running  BISTRO  for  the  12  wheels,  4 layers,  and  44  locations 

of  response  was  2896  sec  on  a CDC  6600,  while  for  AFPAV,  it  was  926  sec 

for  350  elements  and  31  Fourier  terms.  Because  of  the  preprocessor 

associated  with  AFPAV,  the  amount  of  data  preparation  for  both  codes  is 

about  the  same. 

4.2.3  PARAMETER  STUDIES 

Several  parameter  studies  have  been  conducted  to  assess  the  ac- 
curacy of  the  finite  element  predictions  as  well  as  to  provide  a better 
understanding  of  the  effects  of  changes  in  pavement  constituents.  Unfor- 
tunately, this  work  is  usually  an  unreported  part  of  most  studies.  Some 
of  the  work  reported  concerning  understanding  and  accuracy  follows. 
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Figure  22.  C-5A  main  landing  gear  configuration 
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Figure  23.  Response  of  Section  U to  a C-5A  12-wheel  assembly 
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An  extensive  parameter  study  for  both  flexible  and  rigid  linear 
layered  systems  is  reported  in  Reference  19*  This  study  is  primarily 
concerned  with  the  impact  that  changing  the  various  pavement  properties 
has  on  the  predicted  response.  Two  key  points  for  both  rigid  and  flexible 
pavements  are  derived.  First,  soil  stresses  (Figure  23b  is  an  example) 
and  second,  the  displacement  gradients  within  the  deflection  basin  are 
relatively  unaffected  by  significant  changes  in  base,  subbase,  and 
subgrade  moduli.  Figure  2h  demonstrates  that  the  displacement  gradients 
at  the  surface  (i.e.,  the  spacing  between  contour  lines)  remain  relatively 
unchanged  for  the  rigid  two-layer  pavement  systems  shown. 

To  develop  an  accurate  finite  element  prediction  requires  not  only 
a good  mesh  but  also  an  appropriate  soil  island.  This  leads  to  another 
form  of  parameter  study  in  which  the  mesh  and  its  extent  and  displacement 
boundary  conditions  are  examined  with  respect  to  displacement  and  stress- 
strain  convergence.  A convenient  method  for  establishing  convergence  is 
by  comparison  of  finite  element  and  elastic  layer  predictions.  Refer- 
ences 15  and  28  present  such  parameter  studies  with  the  basic  conclusion 
that  a sloping  boundary  of  30  deg  off  the  vertical  with  a free  edge  or 

with  restrained  horizontal  displacement  produces  the  fastest  convergence. 

17 

Unpublished  studies  using  the  WINDAX  Code  indicate  that  finite  element 
displacement:,  within  10  percent  of  those  predicted  by  elastic  layer 
techniques  can  be  achieved  by  using  600  elements  and  soil  island  dimen- 
sions of  100  loading  radii  wide  and  150  radii  deep  for  typical  rigid 
and  flexible  pavements. 

Generally,  stress-strain  predictions  are  unaffected  by  the  size 
and  boundaries  of  the  soil  island  for  any  reasonable  mesh.  Convergent 
displacements  are  considerably  harder  to  obtain  and  are  more  dependent 
on  the  materials  involved.  However,  the  difficulty  associated  with 
displacement  computation  is  not  a serious  liability  for  the  finite 
element  idealization.  First,  to  predict  a linear  displacement  has 
relatively  little  merit  and  as  shown  in  Figures  21  and  23  has  little 
resemblance  to  actual  performance.  Second,  most  useful  displacement 
predictions  are  derived  from  nonlinear  analyses  which  are  much  less 
affected  by  the  soil  island  employed. 
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1*.3  NONLINEAR  FINITE  ELEMENT 


Section  3 presents  a variety  of  material  laws  and  implementation 
algorithms  which  when  coupled  with  the  variety  of  prediction  techniques 
yields  many  alternatives.  However,  in  reviewing  what  is  available  and 
is  being  applied  to  pavements,  there  appear  to  be  only  four  relatively 
widespread  approaches: 

a.  Elastic  layer  analysis  whereby  the  layer  property  is  itera- 

— tively  defined  during  a succession  of  linear  analyses  based 

on  a typical  layer  strain  value.  Such  a method  is  cheap  and 
simple,  and  the  predicted  displacements  are  at  least  nonlinear 
functions  of  load.  However,  as  was  mentioned  earlier,  the 
stress  distribution  and  the  shape  of  the  deflection  basin 
will  remain  relatively  unchanged  as  long  as  the  basic  premise 
of  uniform  layer  properties  is  adhered  to;  this  will  make  any 
precise  usage  of  these  responses  suspect. 

b.  Axisymmetric  finite  element  analysis  using  a one-step  secant 
method  (Figure  15a).  There  are  various  forms  of  this  tech- 
nique available  based  on  the  code  reported  in  Reference  18. 
However,  highly  nonlinear  material  (e.g.,  that  which  cracks 
under  load)  cannot  be  modeled. 

c_.  Axisymmetric  finite  element  analysis  using  the  tangent  method. 
Using  this  technique,  the  response  for  the  repaired  crater 
shown  in  Figure  12  was  computed;  Reference  17  contains  con- 
siderable information  pertaining  to  the  details  of  this  study. 
The  response  of  this  system  where  the  concrete  was  allowed 
and  not  allowed  to  crack  is  shown  in  Figure  25.  In  this 
instance,  the  need  for  a multiple-step  loading,  which  allows 
a cracking  capability,  is  shown  to  be  important  in  ascer- 
taining correlation  between  measured  and  predicted  results. 

£.  Prismatic  finite  element  analysis  using  a one-step  secant 

method  whereby  the  nonlinearity  is  by  element  in  the  X-Y  plane 
and  by  the  layer  along  the  Z-axis  (see  Figure  7).  This  is 
only  a first-order  approximation  to  the  full  nonlinear  analysis 
capability  available  using  the  prismatic  idealization.  Some 
results  are  reported  in  References  21  and  22;  typical  run 
times  are  10-1*0  min  on  a CDC  6600. 

The  axisymmetric  idealizations  offer  the  best  opportunity  to 
study  the  nonlinear  process.  They  are  inexpensive  to  run  (typically 
$3-30);  results  are  relatively  easy  to  interpret;  and  data  preparation 
is  straightforward.  These  codes  provide  an  excellent  environment  in 
which  to  check  the  various  input  parameters  and  evaluate  their  adequacy. 

On  the  other  hand,  the  prismatic  idealization  offers  the  most  realistic 
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procedure  for  modeling  pavements  loaded  by  multiple-wheel  aircraft. 
Also,  the  prismatic  idealization  is  easier  to  "black  box"  for  usage  by 
the  typical  engineer  through  the  usage  of  pre-  and  postprocessing 
software  (Section  5.1*2). 


HORIZONTAL  DISTANCE  FROM  CENTER  OF  LOADED  PLATE,  IN. 


Figure  25.  Deflection  basins  under  30-kip  load  for  repaired 
crater  shown  in  Figure  12 

U.U  THREE-DIMENSIONAL  ANALYSIS 

The  only  relatively  practical  usage  of  three-dimensional  analysis 
known  to  the  authors  is  presented  in  Reference  23.  The  objective  was 
to  make  a rational  determination  of  the  basic  parameters  (e.g.,  pile 
spacing  and  depth)  of  a soil-pile  reinforced  pavement.  Because  of  the 
symmetry  of  the  system  and  the  lack  of  any  other  reasonable  alternatives, 
three-dimensional  analysis  was  used.  Figures  8 and  26  illustrate  some 
of  the  steps  taken  and  the  results  obtained. 
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k . 5 NONLINEAR  KINEMATICS 


Nonlinear  kinematical  considerations  are  required  for  such 
items  as  friction  Joints,1^  mechanically  interacting  multicomponent 
structures  (e.g.,  the  key  Joint  shown  in  Figure  27'3  ),  or  large  deforma- 
tion phenomena  (e.g.,  membrane  surfaces  for  helicopter  landing  zones 

30 

shown  in  Figure  28  ). 


l OCTANT  OF  PILE-REIN FORCEO  b.  MESH  PATTERN  IN  Y-Z  PLANE 

PAVEMENT  MOOEL 


c.  EFFECT  OF  PILE  STIFFNESS  ON  1 VERTICAL  STRESS  AT  INN.  OEPTH 

VERTICAL  SURFACE  OEFLECTION 


Figure  26.  Three-dimensional  analysis  of  a 
pile-reinforced  pavement 
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a.  KEY  JOINT  MODEL 


Figure  21.  Analysis  of  a key  Joint  (from  Reference  38) 
(sheet  1 of  2) 


61* 


Figure  28.  Response  of  membrane  used  as 
surface  (from  Reference  39) 


5 SELECTION  OF  A PREDICTION  TECHNIQUE 


Selection  of  a technique  for  response  prediction  requires  a 
compromise  between  the  Limitations  inherent  in  the  chosen  technique, 
the  cost  of  utilizing  the  technique,  and  its  availability  and  usability. 
Previous  sections  of  this  report  have  alluded  to  most  of  the  limitations. 
The  other  factors  are  discussed  in  this  section  followed  by  some 
recommendations . 

5.1  AVAILABILITY  AND  USABILITY 

Application  of  most  techniques  requires  a computer  code  which  may 
not  be  readily  available;  for  those  that  are,  there  is  still  the  need  to 
adapt  the  code  to  "your"  computer  and  become  familiar  with  its  input  and 
output.  "Teaching  an  old  dog  new  tricks"  was  never  more  applicable  than 
to  pavement  engineers  operating  computer  programs  based  on  various  pre- 
diction theories  which  must  be  coupled  with  performance  and  complicated 
laboratory  and  field  tests. 

The  usability  of  an  analytic  technique  is  ultimately  measured  by 
its  capability  to  aid  in  the  design  of  pavements.  Because  most  tech- 
niques are  integrally  tied  to  computer  hardware  and  software,  it  is  also 
important  that  special  attention  be  focused  on  the  man-machine  interface 
While  the  use  of  computer-based  analytic  procedures  is  still  embryonic 
(especially  finite  element),  two  areas  have  emerged  which  will  play  key 
roles  in  the  future  expansion  of  computer-aided  design.  These  are  per- 
formance models  and  pre-  and  postprocessing  of  analytical  results. 

5.1.1  PERFORMANCE  MODELS 

Although  discussion  of  performance  models  is  beyond  the  scope  of 
this  report,  a brief  comment  on  their  relationship  to  analytic  techniques 
follows.  Computer-based  analytical  techniques  predict  the  displacement, 
stress,  and  strain  responses  of  a pavement  structure  composed  of  various 
idealized  inputs  for  a specific,  one-time  loading.  To  have  meaning  in 
the  design  process,  tnese  responses  must  be  related  to  pav»Tcent  perfor- 
mance, ideally  to  allowable  repetitions  of  loading  prior  t-  * >e  functional 
failure  of  the  pavement  for  a specific  aircraft.  Only  perfo*  ^nce 
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criteria  for  flexible  pavements  have  been  related  to  load  repeti- 
1*0  1*1 

tions,  and  these  are  based  on  relatively  little  data.  Table  2 
presents  those  criteria  known  to  the  authors. 

An  important  "back  door"  approach  to  performance  models  is  usage 
• of  repeated  load  tests  to  derive  material  parameters  which  reflect  the 
worst  possible  condition  of  the  material.  These  parameters,  known  as 
resilient  moduli,  are  input  into  a prediction  technique,  and  voila, 
instant  performance  criteria  are  built  into  the  analysis.  There  are 
limitations  to  such  a procedure:  first,  materials,  usually  in  the 

wearing  course,  that  exhibit  any  brittle  or  fatigue  behavior  (e.g., 
concrete)  cannot  be  treated  this  way;  and  second,  because  soil  materials 
lack  homogeneous  stiffness,  which  aggravates  various  problems  such  as 
asphalt  rutting,  something  more  than  a resilient  modulus  is  required. 
However,  at  this  stage  of  the  game,  such  quantities  as  resilient  modulus 
provide  useful  results  in  an  area  so  lacking  in  other  data. 

5.1.2  PRE-  AND  POSTPROCESSING 
SOFTWARE 

Pre-  and  postprocessing  software  is  an  integral  part  of  computer- 
based  technology  and  to  a certain  extent  is  "the  tail  wagging  the  dog" 
in  that  th%  more  complicated  analysis  techniques,  especially  finite 
element,  are  practically  unusable  without  some  computer  preparation  of 
input  and  computer  generation  of  summary  tables  and  graphical  displays. 
Parenthetically,  techniques  which  are  not  amendable  to  pre-  and  post- 
processing cannot  compete.  Software  packages  can  also  be  developed  which 
provide  combinations  of  techniques.  For  example,  the  BISTRO,  CHEVRON, 
WIND AX,  and  AFPAV  programs  can  be  combined  to  run  using  the  same  input 
and  generating  the  same  type  of  output  where  the  analysis  technique 
is  selected  based  on  the  type  of  problem  to  be  solved  or  the  level  of 
effort  to  be  expended. 

For  preparation  of  finite  element  input,  two  basic  approaches 
exist.  First,  various  levels  of  automation  within  the  finite  element 
context  are  available  (see  References  15,  17,  and  20  for  examples). 

These  primarily  are  schemes  for  generation  of  element  and  nodal  data 
based  on  the  prescription  of  a few  key  elements  and  nodes.  This  type 
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Performance  Criteria  vith  Respect  to  Analytic  Predictions 
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of  automation  probably  requires  an  even  greater  familiarity  with  the 
finite  element  procedure  than  the  relatively  "brute  force"  methods 
associated  with  direct  input  into  the  programs.  For  the  "initiated," 
this  appears  to  provide  the  best  situation. 

The  second  approach  involves  completely  shielding  the  user  from 
the  finite  element  process  (see  References  37  and  42  for  examples). 

Input  is  specified  in  a context  familiar  to  the  design  engineer.  For 
example,  specification  of  a set  of  uniform  layer  thicknesses  is  trans- 
formed by  the  preprocessor  into  a set  of  nodes  and  elements.  This 
approach  offers  the  best  chance  of  incorporating  the  complex  prediction 
techniques  into  the  everyday  design  process.  Paradoxically,  it  is  the 
more  sophisticated  code,  operated  with  fever  a priori  assumptions,  that 
is  more  easily  "black  boxed."  Thus,  making  complicated  procedures  more 
widely  available  is  somewhat  linked  to  making  these  procedures  even 
more  complicated. 

There  are  similar  approaches  relevant  to  postprocessing  of  finite 
element  data  plus  a variety  of  available  hardware  (e.g.,  line  printer, 
drum  plotter,  interactive  and  display  cathode-ray  tube).  Some  graphic 
examples  are  shown  in  Section  4 and  References  15,  17,  and  43. 

5.2  COST  EFFECTIVENESS 

An  important  factor  in  application  of  most  techniques  is  the 
trade-off  between  costs  (of  running  the  technique,  obtaining  the  neces- 
sary input,  and  interpreting  the  output)  and  the  benefits  (derived  from 
a more  refined  or  different  approach).  The  benefits  derived  to  date, 
especially  from  the  finite  element  procedures,  are  largely  educational. 
While  this  has  led  to  a better  understanding  of  how  pavements  perform, 
it  has  not  produced  better  pavements.  To  improve  pavement  designs  with 
any  of  the  more  rigorous  analysis  schemes  will  take  a more  diligent 
and  coordinated  effort  to  establish  the  correlations  between  response 
predictions  and  pavement  performance.  Rudimentary  forms  of  these  cor- 
relations for  flexible  pavement  exist  and  were  presented  in  Section  5.1.1. 

The  cost  of  performing  these  analyses  is  more  clearly  defined. 

For  example,  assuming  all  software  is  in  the  production  mode  and  is 
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being  run  by  experienced  personnel,  a man-week  would  cover  a prismatic 
or  axisymmetric  analysis  excluding  any  detailed  materials  investigation. 
Obtaining  topnotch  laboratory  materials  data  suitable  for  sophisticated 
nonlinear  material  laws  requires  $5,000-$10,000  per  material.  These  are 
the  "bare  bones"  estimates  and  probably  must  be  complemented  with  addi- 
tional dollars  for  training  inexperienced  personnel,  the  necessity  to 
write  and  check  out  new  software,  etc.  Tne  amount  of  these  additional 
dollars  can  easily  overshadow  the  "bare  bones"  figures;  moreover,  most 
design  organizations  are  not  in  positions  to  bring  together  sufficient 
resources  to  make  these  additional  dollars  work  effectively. 

In  sum,  cost  effectiveness  is  not  easily  defined.  Obviously  some 
prediction  technique  is  necessary,  but  whether  today's  methods  or  the 
more  analytic  ones  outlined  in  this  report  are  used  in  future  designs 
cannot  be  settled  by  appealing  to  cost-benefit  arguments  due  to  the  lack 
of  data. 

5.3  RECOMMENDATIONS 

Recommendations  for  selection  of  a prediction  technique  for  typical 
pavement  systems  are  presented  separately  for  linear  layered  pavements 
and  nonlinear  and/or  nonlayered  pavements. 

The  dilemma  of  which  technique  to  use  is  difficult  to  solve. 
Obviously,  a Jointed  runway,  a pile-supported  pavement,  or  a pavement 
crossing  over  a bridge  or  culvert  are  nonlayered  systems.  In  many  cases 
where  experimental  data  have  been  taken,  pavement  response  is  shown  tc 
be  a nonlinear  phenomenon  (e.g..  Figure  23).  However,  the  state-of-the- 
art  does  suggest  that  extensive  use  of  linear  layered  analysis  procedures 
for  design  is  compatible  with  acceptable  practice  and  available  criteria. 
Linear  layered  responses  should  also  be  computed  even  if  only  a prelude 
to  nonlinear  and/or  nonlayered  work.  For  the  near  future,  it  appears  that 
nonlinear  predictions  will  at  least  supplement  linear  ones.  General 
acceptance  of  nonlinear  procedures  seems  dependent  on  establishing  better 
performance  criteria,  making  available  software  (predominately  the  pre- 
and  postprocessing  variety)  better  suited  for  practicing  engineers,  and 
accepting  the  need  for  better  determination  of  material  properties. 
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5.3.1  LINEAR  LAYERED 
PAVEMENTS 

For  a layered  pavement  system  where  linear  responses  are  satis- 
factory, the  elastic  layer  idealization  is  usually  significantly  better 
than  any  other  technique.  Exceptions  would  occur  where  a "hand"  calcula- 
tion is  desirable — possibly  using  the  Westergaard  equations  1&2  - or  when 
a complete  picture  of  response  is  desired  (then  finite  element  would 

usually  be  more  efficient).  Availability  of  the  software  and  hardware 

h-7 

is  a consideration  for  this  method,  but  several  programs  are  rela- 
tively easy  to  obtain.  These  programs  may  be  executed  on  a variety  of 
hardware  for  less  than  $10  for  a "typical"  run.  The  only  actual  handicap 
of  the  elastic  layer  method  is  its  dependence  on  the  accuracy  with 
which  the  layer  properties  are  measured  and  represent  the  field 
conditions. 

The  most  desirable  alternative  to  elastic  layer  calculations  are 
Westergaard  and  axisymmetric  finite  element  idealizations.  As  indicated 
by  Table  1,  the  Westergaard  idealization  produces  results  of  variable 
quality  although  probably  accurate  enough  for  design  calculations  when 
a stiff  layer  tops  relatively  soft  layers.  Its  major  advantages  are 
simplicity  and  usage  of  the  parameter  k . While  k has  some  drawbacks 
associated  with  smearing  all  of  the  soil  properties  into  a single  param- 
eter, it  has  several  important  advantages:  first,  it  is  relatively 

easily  measured  and  is  often  available  from  airport  records;  second,  it 
represents  in  situ  conditions  for  at  least  the  plate  bearing  test  at 
the  time  taken;  and  third,  it  provides  a check  of  elastic  layer  parameters 
(see  Reference  U or  Section  U.l). 

The  axisymmetric  finite  element  idealization,  while  competitive 
in  cost  and  accuracy  to  elastic  layer,  lacks  its  basic  simplicity  of 
operation.  Only  when  appropriate  pre-  «.',nd  postprocessing  software  is 
available,  is  this  method's  use  significantly  more  advantageous  than 
elastic  layer.  Also,  finite  element  usage  generally  requires  added 
expertise  (either  in  software  operation  or  technical  knowledge)  over 
other  methods. 
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5.3.2  NONLINEAR  AND/OR 

NONLAYERED  PAVEMENTS 

It  is  recommended  that  only  those  nonlinear  techniques  which 
allow  variable  material  properties  within  a layer  be  used.  Employing 
this  as  a general  criteria,  the  state-of-the-art  suggests  that  most 
techniques  will  be  based  on  either  the  axisymmetric  (for  a single-wheel 
loading)  or  prismatic  (for  a multiple-wheel  loading)  finite  element 
idealizations  while  occasionally  three-dimensional  idealizations  will 
be  employed  in  peculiar  circumstances  (e.g.  , the  nonlayered  system 
shown  in  Figure  8)  to  predict  primarily  linear  responses.  The  state-of- 
the-art  also  suggests  that  the  material  model  for  soil  be  of  either  the 
secant  or  tangent  variety  (Section  3.1),  and  the  nonlinear  solution 
algorithm  be  of  the  tangent  type  (Section  2.h.2)  and  be  capable  of 
approximating  cracking  and/or  jointed  materials. 
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6 FUTURE  TRENDS  AND  NEEDS 


The  trend  for  pavement  prediction  techniques  is  tovard  more 
sophistication,  accuracy,  complication,  and  automation.  This  is 
generally  true  for  all  forms  of  structural  analysis.  Also,  for  the 
future,  it  appears  that  usage  of  the  finite  element  method  will  in- 
crease significantly.  Some  trends  which  seem  significant  follow: 

a.  New  material  models  for  soils  (e.g.,  Reference  1+1+)  which 
establish  a framework  adequate  for  the  general  situation. 

Also,  modeling  of  other  pavement  materials  will  he  accorded 
more  importance.  Particular  attention  will  be  paid  to  re- 
lating material  mouels  to  performance  criteria. 

b.  New  analytic  procedures,  especially  for  the  finite  element 
method,  will  allow  more  confident  and  efficient  modeling  of 
nonlinear  and/or  nonlayered  pavements.  Use  of  the  so-called 
"global-local"  finite  element  formulation  (References  1+5,  1*6) 
will  help  eliminate  the  uncertainties  of  the  soil  island 
approximation  of  semi- infinite  domains.  This  method  essen- 
tially combines  the  best  features  of  the  finite  element 
method  and  those  of  classical  elasticity.  Finite  elements 
are  used  in  regions  where  loadings  or  high  gradients  of  stress 
occur  while  elasticity  functions — for  example,  those  asso- 
ciated with  elastic  layer  responses,  as  in  equation  3 — are 
used  over  the  whole  idealization.  This  allows  the  finite 
elements  to  be  used  only  where  they  are  needed  and  not  wasted 
merely  to  approximate  the  semi-infinite  nature  of  pavements. 
Another  important  procedure  is  substructuring,  which  provides 
an  efficient  means  for  dealing  with  large  finite  element 
models.  This  allows  portions  of  the  pavement  system  to  be 
analyzed  separately,  ultimately  being  combined  to  produce  a 
unified  result.  For  example,  to  model  joint  intersections 
for  a rigid  pavement,  the  pavement  could  be  idealized  using 
three-dimensional  elements  coupled  along  joint  interfaces, 
supported  by  a separate  three-dimensional  or  prismatic 
idealization  of  the  soil. 

c_.  New  solution  algorithms  for  finite  element  codes  (e.g., 
Reference  1+7)  will  be  more  efficient  and  will  incorporate 
more  advanced  and  accurate  nonlinear  solution  strategies. 

This  will  lead  not  only  to  cost  savings,  but  will  facilitate 
designing  pre-  and  postprocessing  software  that  will  allow 
"black  boxing"  of  the  complicated  features  associated  with 
highly  nonlinear  methods. 

Future  needs  can  be  viewed  from  "technical"  and  "political" 
standpoints.  The  technical  needs  are  relatively  easy  to  define; 
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however , -their 'successful  satisfaction  is  in  part  dependent  on  fulfilling 
the  pdliticai  needs.  Some  important  technical  needs  are  listed  helov: 

a.  The  need  for  standardized,  readily  available,  laboratory 

• materials  testing  procedures  and  equipment- for  establishing 

accurate  constitutive  relations. 

b.  The  need  for  comprehensive  materials  data  available  in  the 
literature  which  is  suitable  for  driving  complex  material 
laws. 

c_.  The  need  for  field  testing  procedures  and  equipment  that  allow 
better  accumulation  of  data  about  in  situ  performance  of 
materials. 

d.  The  need  for  field  data,  especially  measured  strains,  to 
verify  code  predictions  and  material  models. 

e_.  The  need  to  improve  the  usability  of  nonlinear  and/or 

nonlayered  finite  element  techniques.  The  primary  deterrent 
to  widespread  employment  of  these  techniques  is  lack  of  pre- 
and  postprocessing  software  (especially  the  type  that  com- 
pletely shields  the  user  from  the  complex  or  tedious  portion 
of  the  simulation,  see  Section  5.1.2),  and  lack  of  accepted 
performance  criteria  (Section  5.1.1). 

Politically,  the  need  is  to  coordinate  the  pavement  community  in 
both  research  and  field  practices.  Coordination  is  needed  to  deal  with 
problems,  such  as  software  responsibility. 

For  example,  if  a Federal  agency  develops  a design  procedure  based 
on  a particular  computer  code,  how  is  the  code  to  be  made  available  to 
others?  Most  techniques  require  large-size  computers  and  significant 
expense  per  run  and  are  of  sufficient  complexity  to  make  them  relatively 
susceptible  to  misuse.  Thus,  who  is  responsible  for  correctly  installing 
the  code  on  other  computers  and  insuring  its  correct  usage?  Other 
groups  (such  as  the  American  Society  of  Civil  Engineers)  are  studying 
similar  questions.  Some  major  structural  codes  are  widely  distributed 
through  commercial  associations  and  academic  or  government  agencies. 
However,  the  effective  utilization  and  dissemination  of  the  major  re- 
source that  software  has  become  is  still  to  be  achieved. 

Also  needed  are  more  forums  where  the  course  of  pavement  research 
and  practices  may  be  discussed  and  evaluated. 

What  is  aptly  demonstrated  by  the  listing  of  trends  and  needs  is 
the  increasing  dependency  that  design  of  pavements  has  on  a broadening 
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range  of  disciplines.  It  is  not  adequate  just  to  develop  material 
models  without  regard  to  the  availability  of  equipment  for  performing 
the  associated  laboratory  testing,  or  whether  the  law  can  be  incor- 
porated into  an  appropriate  nonlinear  algorithm,  etc.  Those  involved 
in  development  work  should  be  aware  of  where  their  efforts  fit  in  the 
total  effort  related  to  the  prediction  of  pavement  response. 
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